超低轨道卫星应用离子电推进技术方案

低地球轨道大气环境对诸如科学探测和对地观测卫星的阻尼作用十分明显，而且阻尼随太阳和地磁活动以及昼夜、季节交替变化范围宽。为了保证卫星轨道精度或飞行状态满足任务要求，需要利用推进系统对卫星受到的阻尼进行实时或间歇式补偿以实现轨道或飞行状态的保持。针对轨道高度220~268km的无拖曳飞行和轨道维持应用，基于卫星轨道阻尼变化和有效载荷指标要求分析，研究确定了离子电推进技术指标、推力调节方案、系统组成、推力控制方案和在轨应用策略，并对推力调节方案进行了试验验证。结果表明，与无拖曳飞行卫星任务匹配的离子电推进指标为推力调节范围1~20mN，推力分辨率优于12μN，与对地观测卫星轨道维持任务匹配的指标为推力调节范围1～25mN，推力分辨率100μN。研究提出的针对超低轨道卫星应用需求的高精度推力连续调节离子电推进技术方案，具有工程任务针对性和参考价值。     

Contingency target assessment,trajectory design,and analysis for NASA’s NEA scout solar sail mission

The presented study examines contingency target selection and trajectory design for NASA’s Near-Earth Asteroid Scout mission under the assumption of a missed lunar gravity assist. Two previously considered asteroids are selected as potential targets for the given scenario based on favorable orbital characteristics for launch dates ranging from June 27, 2020 through July 26, 2020. Initially, a simplified circular restricted 3-body problem + ideal solar sail model is utilized to survey trajectory options for a month-long launch window. Selected solutions from this data set are then converged in an N-body ephemeris + non-ideal sail model. Results suggest that NEA Scout can still perform asteroid rendezvous mission under the missed lunar gravity assist scenario with new targets, 2019 GF1, 2018 PK21, and 2007 UN12, based on the target launch dates. Further target assessment is carried out for 165 days beyond the current June 27, 2020 launch date.     

多等级区域侦察弹性星座设计方法

针对传统侦察星座目标单一、弹性低的问题，提出了多等级区域侦察弹性星座的设计方法。该方法将星座设计过程按区域等级信息分为多个子星座逐步设计，直到整体星座对所有的区域性能满足设计要求。以区域被划分为3个等级为例，首先对星座设计需求、设计指标及设计步骤进行了分析。其次推导了地面最低分辨率和轨道高度的关系并确定了不同子星座的轨道高度。最后考虑轨道倾角、一箭多星发射、光照和升交点漂移同步约束，构建基础星座、子星座1和子星座2的优化模型。最终设计星座为3层混合星座，共8个轨道面和70颗卫星，星座对各等级区域的最大重访时间分别为10937s，12241s和17437s，弹性指数为2213%，2420%和6361%。结果表明该方法设计的星座可实现对区域覆盖和弹性分级的设计要求，证明了方法的有效性。对比Walker星座设计方法，在同等设计要求下，Walker星座所需卫星数为156颗，多等级区域侦察弹性星座所需卫星数远低于Walker星座，结果进一步证明了该星座设计方法的优越性。     

基于自适应非结构嵌套网格的旋翼流场模拟

针对直升机旋翼CFD仿真的复杂性,提出了改进的适合于格心格式求解器的非结构嵌套网格算法。采用自适应网格技术在旋翼流场仿真的整个过程中进行网格的自适应加密和疏化操作,以更好地捕捉桨尖涡等流动细节。对于频繁的自适应过程中产生的大量重复点和无用点,采用了高效的交替数字树算法(Alternating digital tree,ADT)和标记-删除-移动算法(Mark,delete,move,MDM)进行删除,节约了不必要的存储。针对格心格式的求解器,采用了基于梯度的网格间插值方式,简化了网格间数值传递的复杂性,同时不降低求解器的精度。对CaradonnaTung旋翼悬停算例和HLISHAPE 7A旋翼悬停算例进行了模拟验证,计算值与实验值吻合,表明本文建立的方法具有良好的鲁棒性和有效性。最后,与未采用自适应时求解器对桨尖涡的捕捉效果进行了对比,结果表明本文所采用的方法可以明显地提高求解器对桨尖涡的捕捉。     

基于6σ设计的复合推力高速直升机总体参数多目标优化

针对复合推力高速直升机总体设计阶段总体参数的选择问题,提出一种提高可靠性和鲁棒性的基于6σ设计的改进多目标遗传算法优化方法。采用叶素理论和数值积分的方法分析计算了复合推力高速直升机气动及飞行性能,并以此为基础建立了约束函数和初步目标函数模型;将6σ设计融入改进的多目标遗传算法中,构造最终目标函数;在给定有效载荷设计要求下,对复合推力高速直升机总体参数进行了多目标优化设计。该方法获得了所需的Pareto解,优化后的复合推力高速直升机飞行性能相对原机有了较大改善,算例结果表明该方法有效可行。     

基于概率的运载火箭控制系统设计方法研究

针对传统的控制系统极限偏差设计方法导致运载火箭余量大、总体性能下降的问题,综述了基于概率的控制系统设计方法的国内外研究现状,论述了运载火箭控制系统概率设计的基本思想及设计流程,基于概率密度函数建模理论,建立了通过实际控制指令来控制概率密度函数拟合权值的状态空间,并基于最小二乘方法对状态空间中的参数进行了辨识,进而建立了密度函数成型控制模型,提出了基于最优控制理论的运载火箭控制器设计方法,最后通过仿真验证了理论方法的有效性,为我国运载火箭控制系统的精细化设计提供了理论依据。     

一种基于VPL 的元器件快速建模方法

为了更好地落实航空电子模块的可制造性设计审查工作,结合航空电子模块产品特点和技术要求,建设了基于Valor软件的PCB/PCBA设计与检测自动化分析系统.元器件模型库作为该系统的必要组成部分,其覆盖范围和数据精度对审查结果至关重要.本文将从元器件模型库建设必要性、建设过程和建模工艺要求三个方面进行阐述,并提出一种基于Valor VPL的元器件快速建模方法,应用结果表明,该方法能够高效构建满足系统审查精度的元器件模型,使大规模元器件模型库的建设成为可能.     

空间站燃烧科学实验系统设计

建立空间站燃烧实验系统，可满足未来空间微重力燃烧实验系统需求.通过空间站微重力燃烧实验研究，可拓展空间燃烧学研究.根据所要实现的功能及燃烧实验需求，对中国空间站燃烧柜的燃烧科学实验系统进行了设计和分析.燃烧科学实验系统由8个子系统组成，是一个适合开展气、液、固多种燃料燃烧实验的综合性实验系统.考虑到强度设计要求，在完成方案设计后，对系统进行了有限元分析，并在研制的结构件上进行了力学环境实验.实验与分析结果表明，本文设计的实验系统能够满足环模实验的要求，结构合理可行.      

An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

Multiscale topology optimization using feature-driven method

This paper presents a multiscale design method for simultaneous topology optimization of both macrostructures and microstructures. Geometric features are extended as design primitives at both macro and micro scales and represented by Level Set Functions (LSFs). Parameters related to the locations, sizes, and orientations of macro and micro features are considered as design variables and optimized simultaneously. In the overlapping areas of different macro features, embedded microstructures are optimally figured out as the solution of the corresponding sub-optimization problem. In this study, the eXtended Finite Element Method (XFEM) is implemented for structural and sensitivity analyses with respect to design variables. This method has the advantage of using a fixed grid independent of the topology optimization process. The homogenization procedure is applied to calculate the effective properties of considered microstructures in each macro feature. Numerical examples are presented to illustrate the effectiveness of the proposed method. Results depict that the multiscale design cannot obviously improve structural stiffness compared with a solid-material design under the linear elastic condition.